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Density-functional theory is used to investigate the phase behavior of colloidal binary hard-platelet 
and hard-rod fluids near a single hard wall or confined in a slit pore. The Zwanzig model, in which 
the orientations of the particles of rectangular shape are restricted to three orthogonal orientations, is 
analyzed by numerical minimization of the grand potential functional. The density and orientational 
profiles as well as the surface contributions to the grand potential are determined. The calculations 
exhibit a wall-induced continuous surface transition from uniaxial to biaxial symmetry for the hard- 
rod fluid. Complete wetting of the wall - isotropic liquid interface by a biaxial nematic film for rods 
and a uniaxial nematic film for platelets is found. For the fluids confined by two parallel hard walls 
we determine a first-order capillary nematization transition for large slit widths, which terminates 
in a capillary critical point upon decreasing the slit width. 

PACS numbers: 61.20.-p, 61.30.Gd, 82.70.Dd 



I. INTRODUCTION 



Many complex fluids used in industry or in soft 
condensed-matter laboratories consist of non-spherical 
colloidal particles Q. In particular suspensions of hard 
platelike colloidal particles have recently received ex- 
perimental Ipl, p , ^, ^, 1^ and theoretical attention 
0, !§[ ^, p^, 12], because of the rich phase behav- 
ior and the geophysical and technological implications. It 
has been shown experimentally |5| , theoretically , and 
by simulation [Q, that polydispersity in the size of the 
platelets strongly affects the phase behavior. Whereas 
the theoretical studies have focused on the understand- 
ing of the interactions and the phase behavior of ho- 
mogeneous bulk fluids, experimentally it turns out that 
boundaries such as the walls of the sample cells have 
a pronounced influence on the phase behavior ^ ^. 
Liquid-liquid or wall-liquid interfaces are intrinsic inho- 
mogeneities of the experimental samples which have been 
studied recently |, |, |, ||. 

Here we study inhomogeneous colloidal fluids con- 
sisting of non-spherical particles by examining both bi- 
nary hard-platelet and binary hard-rod fluids within the 
Zwanzig model [|l^ . Platelets or rods are represented by 
square parallelepipeds and the allowed orientations of the 
normal of the particles along their main axis of symme- 
try are restricted to three mutually perpendicular direc- 
tions, rather than a continuous range of orientations in 
space (see Fig. p. Zwanzig's model may be considered as 
a coarse-grained version of the Onsager model which al- 
lows for continuously varying orientations ||l^ . Zwanzig's 
model offers the advantage that the difficult determina- 
tion of inhomogeneous density profiles becomes numeri- 
cally straightforward, allowing one to study various as- 
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pects of inhomogeneous binary hard-platelet and binary 
hard-rod fluids in detail. On the basis of recent experi- 
ence with monodisperse hard-rod fluids |l5|, ^ the 
Zwanzig model is expected to provide a qualitatively cor- 
rect description of the aforementioned colloidal suspen- 
sions by focusing on the entropic properties. In studying 
both binary platelet fluids and binary rod fluids we ad- 
dress the problem of a possible surface transition from 
uniaxial to biaxial symmetry. To the best of our knowl- 
edge properties of binary rod fluids near a hard wall or in 
a slit pore have also not been studied before. Our study 
provides a direct comparison of the structural properties 
and of the behavior of fluids consisting of rodlike and 
platelike particles, respectively. 

In Sec. II we describe the density-functional theory 
and the third-order virial excess free energy functional. 
Section III presents bulk phase diagrams of binary mix- 
tures of thin platelets and binary mixtures of thin rods, 
showing how the density gap at the isotropic-nematic 
transition varies with the mole fraction of the larger par- 
ticles. In Sec. IV we determine the density and orien- 
tational profiles as well as the excess adsorptions of the 
fluids near a hard wall. The calculations exhibit a wall- 
induced surface transition from uniaxial to biaxial sym- 
metry for fluids consisting of rods. Binary hard-rod and 
binary hard-platelet fluids confined by two parallel hard 
walls are investigated in Sec. V. For sufflciently large slit 
widths a first-order capillary nematization transition is 
found. Our results are summarized in Sec. VI. 



II. DENSITY FUNCTIONAL FOR THE 
ZWANZIG MODEL 

We consider a binary mixture of hard rectangular par- 
ticles of size Li X Di X Di {i = 1,2) The number 
density of the centers of mass of the particles of species 
? at a point r is denoted by ^^^^(r). The number den- 
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FIG. 1: The system under consideration consists of a binary 
fluid of thin platelets of surface sizes Di x Di (gray squares) 
and D2 X D2 (black squares) in contact with a planar hard 
wall at z = 0. The figure displays configurations contribut- 
ing to the second order virial term p^\zi) p^J'\z2) [left], and 



the third order virial term px {zi)py {z2)pz (zs) [right] for 
Di = 2D2- For comparison a configuration contributing to 



the second order virial term p^ {zi)py (22) of a binary rod 
fluid is shown. Within this model of only three discrete ori- 
entations particles lying very close to the wall must adopt a 
fully parallel alignment. 



sity of the centers of mass of the particles of species i 
at a point r is denoted by p^p{r). The position of the 
center of mass r is continuous, while the allowed orienta- 
tions of the normal of the particles along their main axis 
of symmetry is restricted to directions (3 = x,y,z. The 
equilibrium density profiles of the mixture under the in- 
fluence of external potentials V^'*-* (r) minimize the grand 
potential functional 



- EE / drpfir) [fcBr(ln[Afp«(r)]-l 



1=1 /3 



pi^ + VPir)] +Fe.[p«(r)] 



(1) 



where are the thermal de Broglie wavelengths and 
Hi are the chemical potentials. Within a third order 
virial approximation the excess free energy functional 
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where fjj^''p^{ri,r2) is the Mayer function. The 
Mayer function equals —1 if the particles over- 
lap and is zero otherwise. With the definition 
•^a /3 = F>i + {Li — Di)Sa.i3, which represents the spatial 
extent in direction a — x,y,z of a particle with orienta- 
tion (3 of the normal, the Mayer function can be written 
explicitly as 



(2) 
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where r^,! is the projection of the position vector ri 
in a direction and 0(r) is the Heaviside step function. 
The density functional theory is completely specified by 
the excess free energy functional and the Mayer func- 
tion. The necessity for including the third order den- 
sity term in Eq. (||), which is not present in the Onsager 
second virial approximation [l^ used in the description 
of thin rods, already follows from recent calculations of 
equilibrium properties of a homogeneous fluid consisting 
of monodisperse thin platelets ^. 

For model systems of hard particles near a structureless 
wall at ^ = 0, apart from a possible surface freezing at 
high densities, nonuniformities of the density occur only 
in the z direction, so that p^p {y) — pp\z). Hence the 
excess free energy functional can be written as 
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where A is the macroscopic surface area in the x — y 
plane. The particular factorization of the Mayer function 
, which results from both particle shape and restricted 
orientations, leads to the relative simplicity of the func- 
tions '^*j'"'/32 (-2^1 J ^2)- Figure |l| displays a second and third 
order virial contribution schematically. 

Third order virial contributions for thin rods (D^ 
Lj, i,j = 1,2) are negligible due to the small intermolec- 
ular interaction between rods [|l3|, |l^ . For thin platelets 
{Li <IC Dj, i,j = 1, 2) the truncation of the virial expan- 
sion after the second order cannot be justified because 
of the nonzero probability of intersection even at small 
thickness Li (see Fig. IT]). 



III. ISOTROPIC AND NEMATIC BULK 
PHASES 

Based on the density functional in Eqs. (|l|), (^, and (j^) 
we study first the homogeneous bulk fluid with V^^\r) — 
in a macroscopic volume V. The equilibrium pro- 

v^y — P/j'') ^n'i the Euler- 
Lagrange equations resulting from the stationarity con- 
ditions dfl[pf]/dpf = for the binary platelet mixture 



files are then constant (^^^(r) 



read ^ ^ [3:,, z ^ j): 
ln(A^p«) = 



M»(fcsT)-i-2Z?f 
~D,Dj{D,+Dj) 



'(J) 

PP2 



.r)2r)4 (i) 0) 



' Ppi 

M 

PP3 



PP2PP3 

PP2 PPz 



(6) 



The Euler-Lagrange equations for a binary mixture of 
thin rectangular rods [D — Di — Dj) are given by ||l8|| : 



ln(Afp«) = p,{kBT)-'-2LjD 



(i) 
Ph 



pV. 



-2L,LjD 



P02 



PPi 



(7) 



The Euler-Lagrange equations (|6|) for the platelet mix- 
ture are independent of the platelet thickness Li because 
the integrals over Mayer functions in Eq. (|^) are indepen- 
dent of Li for thin platelets. For example, the integral 
over the Mayer function of two particles of species i or- 
thogonal to each other is given by 



= y"drl/5^'(rl,0) = 2A(i^ + A)^ 
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In the limit of thin platelets {Li <Si Di) the integral over 

the Mayer function reduces to Vx,'z'' ~ 2Df , which is the 
prefactor of the second term on the right side of Eq. (||). 
In other words, thin platelets have an excluded volume 
(the volume which is denied to a platelet by the condi- 
tion that it must not intersect another platelet), although 
they have a vanishing volume L^-Df (see Fig. H). On the 
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FIG. 2: Relative bulk nematic order parameters si (upper 
curves) and S2 (lower curves) for a binary mixture of thin 
rectangular platelets (solid curves, Di = 2^/^/32) and thin 
rectangular rods (dashed curves, Li = 2L2, Li/D — > 00). 
The number densities of the larger and the smaller particles 
are fixed to pi = P2- Isotropic orientations at low densities are 
characterized by Si = 0, while nematic ordering (0 < Si < 1) 
is preferred for higher densities. 



other hand, the integrals over Mayer functions and the 
excluded volume of thin rods depend on both Li and Di: 
« 2Lf A- Wc have solved Eqs. (|) and (0) numer- 
ically for given chemical potentials. For convenience the 

total number density pf, — J2i ^2/3 P^p ^ ^^'^ number 
densities pi of particles of species i are introduced accord- 
ing to: Pi = px^ +Py^+p^z\ P1+P2 = Pb- The theory has 
been formulated in a way that is completely symmetrical 
with respect to the three coordinate axis. Hence there 
must be a corresponding threefold degeneracy in the re- 
sults. We define the z axis as the preferred coordinate 
axis and consider the relative uniaxial nematic order pa- 



rameters Si = [pz^ — {px^ + py'')/2]/pi. A typical set of 
(si, S2) as a function of pb, with pi = p2, Li = 2L2 for 
rods and Di = 2^/^D2 for platelets is shown in Fig. ||. 
The size ratios of the rods and platelets have been fixed 
such that the second virial coefficients 62 of the equa- 
tion of state of the monodisperse fluids (p2 = 0) in the 
isotropic phase are equal: 



= -Pb [1 + b2Pb + bspl] ksTV , 



(9) 



with 62 = 2LjD/3 for thin rods, and 62 = 2Df/3 for 
thin platelets. For small values of pb the isotropic phase 
(si = S2 = 0) is stable. At a critical density another 
set of solutions, with < Si < 1, appears which repre- 
sents the more favorable nematic phases. The isotropic- 
nematic (IN) transition of the binary platelet mixture 
takes place at a smaller density than for the rod mix- 
ture. We notice that the chemical potentials pi can be 
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FIG. 3: (s) Bulk phase diagram of a fluid consisting of a 
binary mixture of thin platelets (solid curve, 
and a binary mixture of thin rods (dashed curve, Li = 21/2, 
Li/D — > 00) as a function of the chemical potentials /ii and 
fi2 ■ (b) Bulk phase diagram of the same fluids in the density- 
density, i.e., pi — p2 plane (with the same line code as in (a)). 
The straight lines are tielines illustrating isotropic-nematic 
coexistence. In (a) and (b) the dotted lines mark the loca- 
tions of the uniaxial-biaxial transition densities of the binary 
rod fluid in contact with a hard wall. In Figs. |^ and ^ ex- 
cess adsorptions near a hard wall are shown along the two 
thermodynamic paths indicated by arrows. 



written as a function of Si and pi using Eqs. (|^) and (|^). 
At this point it is convenient to introduce the variables 
II* = — fc_B Tin (Af/ci), where Ci = L^Z? for the rod fluid 
and Ci — Df for the platelet fluid. In the following numer- 
ical data are given in terms of /i* and we drop the star in 
order to avoid a clumsy notation. The compositions, den- 
sities, order parameters, and thermodynamic properties 
of the IN coexistence phases are found by solving the co- 
existence conditions: fin = fiiN and pi = where fin, 
fiiN and Pi ~ —ilj/V, pn — —^^n/V are the chemical 
potentials and the pressure of the isotropic and the ne- 
matic phase, respectively. The phase diagrams for binary 
platelet and binary rod mixtures are displayed in Fig. 
The calculations render the concentration in the isotropic 
phase always to be less than in the nematic phase. For 
monodisperse platelet fluids the density gap at the IN 
transition Ap = (p^^-* — = 0.23 is smaller 



than the one for monodisperse rod fluids (Ap = 0.52). 
The different size of the density gap is due to larger in- 
termolecular interactions between platelets as compared 
with those between rods. Moreover, a widening of the IN 
coexistence region is observed at intermediate values of 
the mole fraction of the larger particles. With increasing 
the size ratios of the particles (Li = 3L2, Di = 8^/^1)2) 
the calculations exhibit a greater degree of fractionation 
between the two coexisting phases. The smaller particles 
are preferentially in the isotropic phase. The results for 
rods are in agreement with earlier calculations eM. 



IV. HARD-ROD AND HARD-PLATELET 
FLUIDS NEAR A HARD WALL 

The density and orientational profiles of both compo- 
nents of binary hard-rod and binary hard-platelet mix- 
tures close to a planar hard wall are obtained by a nu- 
merical minimization of the grand potential functional 
(Q) with the excess free energy functional (^) . The results 
are conveniently expressed in terms of the orientationally 
averaged number density profiles 



p,(z)=p«(z) + p«(z)+p«(z). 
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position-dependent nematic order parameters 
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and position-dependent biaxial order parameters 



(11) 
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At small distances from the wall the value of the ne- 
matic order parameter reflects the geometric constraints. 
A platelet [rod] lying very closely to the wall must adopt a 
fully parallel alignment (see Fig. |l]), so that the nematic 
order parameters reach their limiting values Si(0) = 1 
[si{0) = —1/2] there, whereas the isotropic orientation 
Si{z) = is attained at large distances from the wall. 
Orientational profiles of biaxial symmetry are described 
by Si{z) 0, qi{z) ^ 0, where a positive or negative sign 
of qi{z) signals a spontaneous preferential alignment of 
the normals parallel to the x-axis or y-axis, respectively. 



A. Monodisperse fluids 

The phase behavior of monodisperse hard-rod fluids 
(^2(2) — 0) near a structureless wall has been investi- 
gated using the Zwanzig model and a wall-induced tran- 
sition from uniaxial to biaxial symmetr y up on increas- 
ing the bulk density has been found |]l5[ |l^. Moreover, 
we have recently investigated both monodisperse hard- 
platelet fiuids and monodisperse hard-rod fluids near a 
structureless wall at low bulk densities using a model that 
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FIG. 4: Orientationally averaged number density profiles 
pi{z), nematic order parameters si(z), and biaxial order pa- 
rameters qi[z) for monodisperse thin platelets (solid curves) 
and thin rods (dashed curves, D/Li ^ oo) in contact with a 
planar hard wall at « = 0. Positive [negative] values of the 
nematic order parameters indicate that the platelets [rods] 
are preferentially aligned parallel to the wall. Biaxial sym- 
metry of rods near the wall is characterized by qi{z) 7^ 0. 
Here qi{z) 7^ for z > Li/2. The bulk density is fixed to 
piDi = 1.13 for the platelets and piL\D = 1.13 for the rods. 



profile and of the order parameters: 

ln[A?p(i)(z)]-ln[A?pW(z)] 
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for the hard-rod fluid and 
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for the hard-platelet fluid. One easily finds that the 
uniaxial distribution, with qi{z) 0, is a solution of 
Eqs. ([ij) and ( |l5| ) for any density pi(z) and nematic 
order profile si(z). Biaxial distributions qi{z) ^ for 

the rod fluid are possible if pi{z) > pf^^\z), where the 

uniaxial-biaxial (UB) transition density p'f^^'^ {z) follows 
from a low-gi expansion of Ai(2:) = 3gi(z)/[l — si(z)] + 
0{ql{z))iov si(z)^l: 



{UB), ^ 



2[l ~ si{z)]LlD 
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Since si{z) — ~\ /2 {Q < z < Li/2) is the minimum value 

of si{z) for rods and because p^^^\z) decreases with 
decreasing si(z), it follows that local biaxiality starts to 
develop if pi{z) = 1 in the interval < z < Li/2 (see 
Fig. |). For platelets the value si(^) = 1 (0 < z < Di/2) 
is determined by the geometric constraint. Hence qi{z) — 
is the only solution of Eqs. ([l^) and (|l5|) close to the 
wall. 



allows for continuous orientations of the particles [ p2| . 
Here we extend these previous calculations and present 
in Fig. I the calculated order parameters of monodisperse 
hard-platelet fluids and monodisperse hard-rod fluids for 
high bulk densities. The sharp cusps at z — Di/2 for 
platelets, and at ^ = Li/2 for rods, reflect the disconti- 
nuities of p''x\z) and p^\z) for platelets, and of p^z\z) 
for rods, which determine the value of the nematic order 
parameters close to the wall. The most noteworthy fea- 
ture is that no biaxiality is found for the platelet fluid 
{qi{z) ~ 0) while the loss of translational invariance due 
to the hard wall breaks the uniaxial symmetry of the 
rod fluid {qi{z) 7^ 0). In order to study the possible on- 
set of biaxiality, we rewrite the Euler-Lagrange equations 
in terms of the orientationally averaged number density 



B. Binary fluids 

The binary fluids considered in this subsection consist 
either of two types of thin rods (Li = 2L2) or of two 
types of thin platelets {Di = 2'^/^D2)- We focus on the 
numerically determined excess adsorptions defined as 



r, = 



dz [pi[z) - Pi 



(17) 



where pi — pi{z 00). Figures |^ and ^ display F^ 
for binary rod fluids and binary platelet fluids, respec- 
tively. For a fixed bulk density of the small particles p2 
and for small bulk densities of the large particles pi , the 
excess adsorption of the small particles increases upon 



increasing pi. The reason for this is that the increas- 
ing number of large particles lying close to the wall en- 
forces the orientational ordering of the small particles, 
leading to an enrichment of small particles near the wall 
because of reduced intcrmolccular interactions between 
the latter as compared to the isotropic bulk fluid. Due 
to the same mechanism, Fi increases upon increasing p2 
for constant bulk densities of the large particles pi . The 
excess adsorption of the large particles exhibits a change 
of sign and sharp increase with increasing pi while a net 
depletion of the small particles is found for small bulk 
densities p2- The calculation renders Fi to diverge log- 
arithmically as pi — > p^^\ where p^^^ is the bulk den- 
sity of the large particles in the isotropic phase at the 
IN transition (see lower curves in the upper figure of 

Fig. ^). Near p^^^ the excess coverage can be fitted by 
Fi = Ai - A2ln(L2£)[pW _ p^]) for the rod fluid and 

Fi = Si - B2 \yi{DI[p^^^ - pi]) for the platelet fluid, with 
fit parameters Ai, A2 and Bi, B2, while the excess cov- 
erage of the small particles remains finite and attains 



the critical value F, 



.(c) 



consistent with complete wetting of the wall 



via a square-root cusp singularity 

The logarithmic increase of Fi is 
isotropic 



fluid interface by a nematic fllm at p\ . Complete wet- 
ting is confirmed explicitly by the vanishing of the contact 
angle (see Sec. IV) . One observes that the UB transition 
of the rods marks the onset of a pronounced variation of 
Fi, while Fi increases more smoothly upon increasing pi 
for the binary platelet fluid, due to the absence of biax- 
iality. In order to study the possible onset of biaxiality 
for the binary rod fluid, we rewrite the Euler-Lagrange 
equations in terms of the density profiles and order pa- 
rameters: 



In 



1 - si{z) + |gi(z) 



l-.si(z)-|gi(2) 



L2 



In 



1 - S2(z) + §92(2) 



1-S2{Z) 

2D[Llpi{z)qi(z) 



192(2). 
2Llp2{z)q2{z)] 



(18) 



For any Si{z) and Pi{z) there is a trivial uniaxial so- 
lution qi{z) = to Eq. (p^). Nontrivial biaxial solu- 
tions qi{z) ^ exist at sufficiently high local densities 
Pi{z) > p[^^\z) with the UB transition densities follow- 
ing from an expansion of the logarithms in Eq. (|l8|): 
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2[l-si(z)] 

>iz)LiD 
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In the case of a monodisperse hard-rod fluid Eq. ( pj| ) 
reduces to Eq. (|l^). Since si{z) — S2{z) — —1/2 close to 
the wall, local biaxiality sets in if pi{z) Li^/ L'\p2{z) = 1 
near the wall. We notice that, independent of the size 
ratios of the rods, biaxiality of one species of the fluid is 
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FIG. 5: The excess adsorptions Fi and r2 of large ((a)) and 
small ((b)) rods, respectively, (Li = 2L2, Li/D 00) near 
a hard wall as a function of the bulk density p\ of the large 
rods for two values of the bulk density of the small rods: 
P2LID = 2 (solid lines), P2LID = 0.5 (dashed lines). The 
corresponding thermodynamic paths are indicated in Fig. ^ 
(a). The diamonds and the vertical lines mark the location 

of the uniaxial-biaxial transition densities Pi*^^' and the den- 
sities of the isotropic bulk phase at isotropic-nematic coex- 
istence Pi^\ respectively. For comparison, p^^^' and p2 are 
marked by the dotted line and arrows in the bulk phase dia- 
gram displayed in Fig. |^ (a). Fi diverges logarithmically as 
pi — > Pi^', while r2 attains finite values (*) via square- root 
cusp singularities. 



always accompanied by biaxiality of the other species, as 
expected on geometrical grounds. 



V. BINARY HARD-ROD AND 
HARD-PLATELET FLUIDS CONFINED BY TWO 
PARALLEL HARD WALLS 

The results of the previous section show that a hard 
wall favors planar nematic order (the main body of the 
particles is oriented parallel to the wall) over isotropic or- 
der upon increasing the particle densities. We now con- 
sider binary hard-rod and hard-platelet fluids confined 
by two parallel hard walls at z = and z = H and in- 
vestigate a possible capillary condensation of a nematic 
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FIG. 6: The excess adsorptions Fi and F2 of large ((a)) and 
small ((b)) platelets, respectively, {Di = 2^^^D2) near a hard 
wall as a function of the bulk density pi of the large platelets 
for two values of the bulk density of the small platelets: 
P2-D1 = 2 (solid lines), p2-D? = 0.5 (dashed lines). The cor- 
responding thermodynamic paths are indicated in Fig. ^ (a) . 
The vertical lines mark the locations of the densities of the 
isotropic phase at isotropic-nematic coexistence p['\ Fi di- 
verges logarithmically as pi Pi \ while F2 attains finite 
values (*) via square-root cusp singularities. 



phase. Particularly, we calculate the surface contribu- 
tions defined via 



where A is the area of a single surface, Wf, is the bulk 
grand-canonical potential density, and V is defined as 
the volume of the container with its surface given by the 
position of the rim of the particles at closest approach so 
that V = AH. 7u;7 is the wall - isotropic liquid surface 
tension in the absence of the second wall and ^(i/) is 
the finite size contribution. We restrict our attention to 
chemical potentials fii smaller than the chemical poten- 
tials /J-i^^^ at bulk isotropic-nematic coexistence. 



A. Surface tensions and wetting at a single hard 
wall 

Figure ^ displays the surface tension 7^,/ as a func- 
tion of the chemical potential of the larger particles. The 
steric interaction between the particles, which is more 
pronounced for the platelets, increases the surface ten- 
sion with increasing chemical potential. On the other 
hand, nematic ordering of the particles, induced by the 
walls, leads to a decrease of the surface tension for large 
chemical potentials. For fixed small chemical potential 
fi2 of the small particles a maximum of 7^7 as a func- 
tion of the chemical potential fj,i of the large particles is 
observed. 

For large negative chemical potentials, i.e., in the limit 
of noninteracting particles, the wall - isotropic liquid sur- 
face tension is given by 



jwiLiD _ 1 



for the binary rod fluid, and 
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(22) 



for the binary platelet fluid. The prefactor 1/2 in 
Eq. (^l|) reflects the fact, that the orientationally aver- 
aged excluded volume due to the wall is smaller for rods 
than for platelets. 

For the binary rod mixture 7^,7 is a non-monotonic 
function of fii close the uniaxial-biaxial transition. The 
local minimum of 7^,7 displayed in the inset of Fig. |^ 
disappears upon increasing or decreasing fi2, i.e., in the 
limit of monodisperse fiuids. 

We have confirmed that complete wetting of the wall - 
isotropic liquid interface by a nematic film occurs along 
the whole isotropic-nematic coexistence by observing a 
vanishing contact angle d: 



cos-i? 



Here 7ui/(Mi^^'*) is the wall - isotropic liquid surface ten- 
sion, 7M)Af (Mj-^^'') is the wall-nematic liquid surface ten- 
sion, and 77 AT is the isotropic-nematic interfacial tension. 
All tensions are taken at isotropic-nematic two-phase co- 
existence. The chemical potentials at the IN transitions 
are denoted by ^-^^^ (see Fig. || (b)). 



lIN 



(23) 



B. Film geometry 

The results for the finite size contribution ui{H) are 
shown in Fig. ^. For sufficiently large slit widths H the 
slope of ll!{H) as a function of the chemical potential 
of the large particles changes discontinuously at a criti- 
cal value, signalling a first-order capillary nematization 




FIG. 7: The wall - isotropic liquid surface tension 7^7 of a 
binary hard-rod fluid (Li — 2L2, L\/D 06), (a), and a bi- 
nary hard-platelet fluid (Di = 2^/^1)2), (b), respectively. The 
chemical potential of the small particles is kept flxed for each 
curve and increases from bottom to top: fi2/{kBT) = —1, 0, 1. 
Bulk isotropic-nematic coexistence occurs at /ii = ^\^'^\jj,2), 
which corresponds to the final points of the curves at the right 
(see Fig. |^ (b)). The inset displays 7^,/ of the binary hard- 
rod fluid for 112 /(kBT) — 1 with increased resolution. The 
diamonds mark the uniaxial-biaxial transitions. 



transition. We emphasize that lij{H) is a smooth and 
monotonic function of /ii for fixed H close to the uniaxial- 
biaxial transition whereas as discussed above 7^,7 exhibits 
a local minimum for binary rod mixtures close to the 
uniaxial-biaxial transition. As function of H the finite 
size contribution uj{H) corresponds to the solvation free 
energy for the immersed two plates acting as the confin- 
ing walls for the fluid. According to the insets in Fig. ||, 
iij{H) exhibits a minimum at iJ = (by construction, 
Lo{0) = — 27^,7 and Ld{oo) = 0) so that the solvation 
force —du){H)/dH is attractive. For a discussion on the 
repercussions of this solvation force in the context of col- 
loidal stability we refer to, e.g., Ref . J9|| . Whereas for 
the present systems > (see Fig. |3)7for simple fluids 
the surface tension 7^; is typically negative whenever the 
wall-liquid attractions dominate over liquid-liquid attrac- 
tions and even for hard walls ^wi < at densities well 
above liquid-vapor coexistence for which drying does not 
occur. In those situations w(0) > and the solvation 



FIG. 8: The flnite size contribution ijj{H) of the grand po- 
tential functional [see Eq. (po|)] of a binary hard-rod fluid 
(Li — 21/2, L\/D 00), (a), and a binary hard-platelet 
fluid (Di — 2^^'^D2), (b), respectively, confined in a slit of 
width H , and in contact with an isotropic bulk reservoir at 
chemical potential ^1. The chemical potential of the small 
particles is kept fixed at /i2/(fcsr) = —1. The width of the 
slit increases from top to bottom: H/L\ = 1.5, 1.75, 2.5, 3 in 
(a) and H / Di = 3, 3.5, 5, 6 in (b). Bulk isotropic-nematic co- 
existence occurs at ^\ — fi['' '^\fi2) , which are the maximum 
/xi -values displayed (see Fig. ^ (b)). The insets display the sol- 
vation free energy oj{H) as function oi H at /7i/(fcsT) = 0.64. 



potential decreases as H increases. 

The occurrence of the capillary condensed nematic 
phase can be inferred more directly from the density and 
nematic order profiles displayed in Fig. ^. The density 
profile of the capillary condensed nematic phase is char- 
acterized by a nematic phase throughout the slit, whereas 
the density profile of the coexisting phase decays toward 
an isotropic phase in the middle of the slit. As expected 
from the presence of purely repulsive walls, the total lo- 
cal midplane density pi{H/2) + p2{H/2) of the capillary 
condensed nematic phase is slightly smaller than the co- 
existing nematic bulk density p^^^ -I- ■ For the binary 
rod fluid both the capillary condensed nematic phase and 
the nematic film in the coexisting isotropic phase are bi- 
axially symmetric. The location of the uniaxial-biaxial 
transition has practically not been altered by the con- 
finement. 



It is apparent from Fig. ^ that the interfacial profile 
in the isotropic phase is larger for the platelet fluid than 
for the rod fluid because of the relatively smaller inter- 
molecular interactions between rods as compared with 
those between platelets. We observe coexistence between 
the isotropic and the capillary condensed ncmatic phase 
provided H > iJc(/^i,M2)- For sufficiently narrow slits 
(H < iJc(/^i,M2)) a sharp capillary nematization tran- 
sition no longer occurs and is replaced by a steep but 
continuous filling. Hence the capillary nematization tran- 
sition ends in a capillary critical point at a critical wall 
separation iJc(/xi, ^2)- 

We have determined the capillary nematization tran- 
sition for various slit widths and the phase diagrams for 
binary rod and binary platelet mixtures constructed as 
function of the chemical potentials and the slit widths are 
shown in Figs. |l^ (a) and |ll| (a). Upon decreasing the 
slit width (H > Hc{^i, ^12)) the capillary nematization 
transition is shifted to smaller chemical potentials rem- 
iniscent of the shift of the capillary condensation tran- 
sition in confined simple liquids. Upon increasing the 
chemical potential of the small particles the critical 
wall separation decreases because the small particles are 
preferentially in the isotropic phase in the middle of the 
slit. This leads to a depletion of the large particles in the 
middle of the slit and prevents a continuous filling of the 
slit even for rather small slit widths. 

Figures |l^ (b) and |ll| (b) display an alternative rep- 
resentation of the capillary phase diagrams in terms of 
average number densities defined as 
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H 



dz pi{z) 



(24) 



Here the curves with low values of {pi) represent the 
coexisting isotropic phase and the curves with high val- 
ues of {pi) the capillary condensed nematic phase. Upon 
approaching the critical point (H < Hc{{pi) , (P2))) the 
difference of the average densities in the coexisting phases 
becomes smaller and vanishes at the critical point. The 
shape of the coexistence curves in Figs. |l^ (b) and |ll| (b) 
reflects the fact that the difference in {pi) between the 
two phases is used as order parameter. Our calculations 
lead to mean-field order parameter exponents; hence the 
dash-dotted lines in Figs. |l^ (b) and |l^ (b) do not inter- 
sect under an acute angle. The actual rounding becomes 
visible only at higher resolutions. 



VI. SUMMARY 

We have applied a density-functional theory to fluids 
consisting of binary hard-platelets and hard-rods near 
a single hard wall or confined in a slit pore of width 
H . The particles are square parallelepipeds with orien- 
tations restricted to three mutually perpendicular direc- 
tions (Fig. 1^). The rectangular shape of the particles and 
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FIG. 9: Coexisting isotropic and capillary condensed nematic 
profiles for a binary hard-rod fluid (L\ = 2L2, L\/D 00), 
(a), and a binary hard-platelet fluid (Di — 2^^^D2), (b), con- 
fined in a slit of width H = 6L1 and H = 6D1, respectively. 
The isotropic profiles of the large particles (dashed lines) are 
characterized by si{z) = in the central region, while neg- 
ative [positive] values of the nematic order parameters si{z) 
close to the confining walls at 2: = and z = 6L1 [z = 6-Di] in- 
dicate that the rods [platelets] are aligned parallel to the walls. 
The capillary condensed nematic phase (solid lines) is charac- 
terized by strong orientational ordering {si{z) 7^ 0) through- 
out the slit. The chemical potential of the small particles is 
P2/{kBT) = -1.0 so that in the bulk p['^ LID = 0.44 and 
p['^LlD = 1.14 in (a) and pi'^ Df = 0.44 and p['^ Df = 1.02 
in (b). 



their restricted orientations allow one, within the frame- 
work of a third-order virial approximation of the excess 
free energy functional, to determine numerically the den- 
sity proflles, orientational profiles, surface and finite-size 
contributions to the grand potential, and phase diagrams 
with the following main results. 

(1) Figures || and || demonstrate that the bulk 
isotropic-nematic transition of the binary platelet mix- 
ture (size ratio: Di = 2^/^D2, where Di x Di is the 
surface size of the thin platelets of species i) occurs at 
a smaller density than for the rod mixture (size ratio: 
Li — 2L2, where Li is the length of the thin rods of 
species i). Moreover, the density gap at the isotropic- 
nematic transition is smaller for the rod fluid than for 




FIG. 10: (a) Phase diagram of a confined binary hard-rod 
fluid (Li — 21/2, Li/D — > cx) as a function of the chemi- 
cal potentials fii, fi2, and wall separation H. The solid line 
represents the bulk phase diagram (see Fig. ^ (a)), while the 
dashed and the dash-dotted lines correspond to wall separa- 
tions H = 61/1 and H — 2L\. For fixed H < Hc{^jli, 112) 
the capillary condensation transition is first order and termi- 
nates at a critical point. These critical points for various slit 
widths form the dotted line. One example of such a capillary 
critical point is indicated by the solid circle. There the corre- 
sponding line of first-order capillary transitions for H = 2Li 
(dash-dotted line) ends, (b) Phase diagram of the same fluid 
plotted (with the same line code as in (a)) as function of the 
average number densities (pi) and {p2) in the slit. In between 
corresponding lines there is two-phase coexistence between 
isotropic and capillary condensed nematic phases. For a small 
slit width H — 2Li (dash-dotted fine) the branches of the 
coexisting capillary condensed nematic phase and isotropic 
phase end at a critical point denoted by the solid circle. 



the platelet fluid because of larger intermolecular inter- 
actions between platelets as compared with those of rods. 

(2) Platelets [rods] lying very closely to a planar hard 
wall must adopt a fully parallel alignment (see Fig. 0) , so 
that the nematic order parameters (Eq. (O)) reach their 
limiting values Si(0) = 1 [s^(0) = -1/2] there (Fig. |). 
The biaxial order profiles (Eq. (|l2[)) for the inhomoge- 
neous fluids at densities slightly below the bulk isotropic- 
nematic transition densities demonstrate biaxial symme- 
try {qi{0) 7^ 0) of rods and uniaxial symmetry {qi{0) = 0) 
of platelets near the wall (Fig. 



FIG. 11: Phase diagrams of a confined binary hard-platelet 
fluid {Di — 2^^^ 02). The solid lines represent the bulk phase 
diagram (see Fig. pi), while the dashed and the dash-dotted 
lines correspond to wall separations H = 6D1 and H — 2Di. 
The critical points for various slit widths form the dotted line 
in (a). 



(3) The excess adsorption Fi of the large particles di- 
verges logarithmically as pi — s- p[^^ , where p\^^ is the bulk 
density of the large particles in the isotropic phase at the 
isotropic- nematic transition (Figs. ^ and |6[). Complete 
wetting of the wall - isotropic liquid interface by a ne- 
matic film is confirmed by observing a vanishing contact 
angle. 

(4) The steric interaction between the particles in- 
creases the wall - isotropic liquid surface tension 7^,/ 
with increasing chemical potential, whereas nematic or- 
dering of the particles, induced by the walls, leads to a 
decrease of the surface tension for large chemical poten- 
tials (Fig. 0). 

(5) For sufficiently large slit widths the slope of the fi- 
nite size contribution to the grand potential as a function 
of the chemical potential of the large particles changes 
discontinuously at a critical value, signaling a first-order 
capillary nematization transition (Fig. H). The density 
profile of the capillary condensed nematic phase is char- 
acterized by a nematic phase throughout the slit, whereas 
the density profile of the coexisting phase decays toward 
an isotropic phase in the middle of the slit (Fig. |^). 
The isotropic-nematic interfacial profile is broader for the 
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platelet fluid than for the rod fluid. 

(6) Coexistence between the isotropic and the capillary 
condensed nematic phase is observed provided the slit 
width H is sufficiently large: H > Hc{iJ.i, 1^2)- For suf- 



ficiently narrow slits {H < Hdni, H2)) a sharp capillary 
nematization transition no longer occurs and is replaced 
by a steep but continuous filling (Figs. ^ and [Tl] ). 
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